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The  Extinction  Probability  in  a 


Critical  Branching  Process 
by  Howard  J.  Weiner* 
University  of  California  at  Davis 
and  Stanford  University 


I . Introduction. 

Let  Z(t)  denote  the  number  of  cells  alive  at  t in  a critical  age- 
dependent  Bellman-Harris  branching  process  with  cell  lifetime  distribution 
G(t),  G(0+)  = 0,  non-lattice  and  assume 


(1.1) 

(1.2) 


t (l-G(t))  -*  0 as  t -* 

CO 

0 < \s.  s J tdG(t)  < ® . 


0 


Let  the  offspring  generating  function  be  denoted,  for  0 < s < 1, 

CO 

k 


(1.3) 


and 


(1.4) 


and 


h (s)  = £ p s 

k=0 


h (1)  = 1 — £ kp 

k=l  * 


(criticality) 


(1.5) 

See  [1]  Chapter  4 for  details. 


0 < a2  s h"(l)  = £ k(k-l)P  < 

k=2  k 
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It  is  well-known  that 


(1.6) 


lim  tP[Z(t)  > 0]  = ~ h . 

p — » CO 


Various  proofs  of  (1.6)  and  the  corresponding  result  for  a critical 
Galton-Watson  or  discrete  time  process  have  appeared.  See  [1]  Chapters  1,  4 
for  comments  and  references.  For  example,  the  proof  of  (1.6)  in  [2]  uses  the 
renewal  theorem.  The  proof  in  [1] , Chapter  4,  gives  the  asymptotic  form  for 
a critical  generating  function,  from  which  (1.6)  is  a special  case,  and 
relates  the  generating  function  to  that  of  a critical  Galton-Watson  process, 
for  which  results  are  obtained.  The  proof  given  in  this  note  is  elementary 
and  self-contained. 


II . Comparisons  and  Iterations 

Definition.  For  0<s<  1,  t>0. 


(2.1)  F(s,t)  = 2 P[Z(t)=k]sk. 

k=0 

It  is  well-known  ([1],  Chapter  4),  and  follows  by  the  law  of  total 
probability  that 

(2.2)  F (s , t)  = s (l-G(t) ) + J h (F (s , t-u) )dG(u)  . 

0 

Using  the  fact  that 


(2.3) 


P(t)  3 P[Z(t)  >0]  = 1-F (0,t)  , 


then  from  (2.2)  it  follows  that,  by  a Taylor  expansion  on  h(s)  about  s=l , 


3 


(2.4) 


and 


-/ 


l-P(t)  = j h(l-P(t-u) )dG(u) 
' 0 


t 2 

(2.5)  l-P(t)  = J {l-P(t-u)  + ~ P2(t-u)  + o (P2 (t-u))}dG(u) 


Rewriting  (2.5), 

r a2  rfc  2 

(2.6)  P(t)  = l-G(t)  +J  P(t-u)dG(u)  -j  J P (t-u)dG(u)  +f(t) 

~ 0 


where  f(t)  denotes  the  remainder. 


Lemma  1 ([!])  P(t)|0. 


Proof.  A simpler  and  elementary  proof  will  be  given  here.  Clearly 


(2.7) 


P(t)  i C > 0. 


Assume  C>0.  Split  the  integral  on  the  right  side  of  (2.5)  into  J = 
pfc/2  rt  _2  o 

+ J . For  large  t,  the  integral  J is  o(t  ).  Then  it  is  clear  that 
0 t/2  t/2 


(2.8) 


1-C  - 1-C+-  C+o(l), 


which  is  a contradiction  of  C>0.  This  proves  the  lemma. 


(2.9) 


Let  G^n^  (t)  = nt^'  iterate  of  G evaluated  at  t. 


Define  the  iterative  sequences 


(2.10) 


Un+l(t)  = I h(Un(t-u)dG(u) 


U (t)  s 1. 
o 


4 


(2.11) 


1 ,1 = f hd  (t-u))dG(u) 
n+1  J _ n 


A 


i (t>  » 

o 


, t < T 


V. 


1 - -3  t > T 


for  T > b. 


Lemma  2 Under  assumptions  (1.1)  - (1.5),  for  n > 0,  and  all  t > T, 

(2.12)  0 < U (t)  - F(0,t)  < G(n)(t) 

(2.13)  0 < U (t)  - I (t)  < G(n)(t) 

— n n — 


(2.14) 


In(t)-Io(t) j < k(t). 


where,  for  t -*  «>, 


(2.15) 


tk(t)  - 0. 


Proof.  Eq. 
Assume 


(212)  holds  for  n = 
(2.12)  for  n.  Then, 


0 . 

omitting  arguments. 


(2.16) 


0 < Un+1-F  = J h(Un)-h(F)dG  < J (U^-F)dG 


<J  G(n)(t-u)dG  = G(n+1)(t) 


where  the  left  inequality  follows  from  the  induction  hypothesis  and  the 
monotonicity  of  h,  and  the  right  inequalities  from  the  mean  value  theorem,  the 

t 

fact  that  h (1)  =1,  and  the  induction  hypothesis. 

To  show  (2.13),  observe  that  for  t >T 


0 < U -I  < - < 1 
— o o — t 


G(o)(t)  . 


(2.17) 


Assume  (2.13)  for  n by  induction..  Then,  arguing  as  in  (2.16) 


(2.18)  0 < U -I  = f (h(U  )-h(I  ))dG  < f(U  -I  )dG 

— n+1  n+1  J n n n — J _ n n 


pt  .(n) 


< j GyilJ  (t-u)dG(u)  = G^n+1^(t)  . 

J0 


To  show  (2.14)  write  I (t)  as 


pt  / 2 pc 

I. (t)  — J + h (I  (t-u) )dG(u) 

1 - t/2  0 


(2.19) 


0 


-2 


Note  that  the  second  integral  (J  ) in  (2.19)  is  o (t  "”)  by  (1.1). 

pt/2  t/2 

The  first  integral  ( ) may  be  written  by  a Taylor  expansion  of  h(s)  about 


s = 1 as  (t  » T) 

(2.20) 


0 


2 Pt/2  h 

1 (t)  = o(t  ) + J h(l-~—  )dG(u) 

i 0 


o pt/2  1 K 9 T * 

o(t  ) + J {h(l)  - (~^)h'(l)  + 7(7^7)  h (1) 


2 t-u 


+ o(t~2)}  dG(u) 

! 

Using  h(l)  = h (1)  = 1,  and  the  expansion,  0 < u < t/2. 


(2.21) 


— = - (1  +— +o(— ) ) 
-t-u  t t vt 


in  the  right  side  of  (2.20)  yields 


t / 2 

(2.22)  L(t)=o(t'2)+G(t)-;f  (1  +7  + o(J)  ) dG(u) 

0 

+ ^2  J 1 ^ +o(t)  * dG(u) 

2t2  0 

as  again  by  (1.1),  an  integration  by  parts  yields 

CD 

(2.23)  J udG(u)  = o(t  d) . 
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It  follows  that 


(2.24) 


IL(t)  = IQ(t)  + f(t) 


where 


(2.25) 


0 < Jf  (t)  j < K < co 


and  as  t -»  00 , 


(2.26) 


t2f(t)  - 0. 


Then  one  may  write 


(2.27) 


I2(t) 


h(In(t-u)+f (t-u))dG(u) 

l)  ~ U 


> 


[h(I0(t-u))+h' (IQ(t-u))f (t-u)}dG(u) . 


Note  that  for  t » T, 


(2.28) 


T 

(T)  = J h(I0(T-u))dG(u) 


and 

(2.29)  h(l-e)G(T)  < I1(T)  < h(l  - |)G(T), 


for  some  a > 0S  e > 0, 


(2.30)  I2(t)  < J Ch(I0(t-u)  + h ' (max{h (1  - |)G(T),1  - ^})f 


t-u) }dG(u) 


for  some  a > 0,  y > 0. 


7 


For  all  t » T,  (2.27)  - (2.30)  yield  that 


ft 


(2.31) 


{ I2(t)-I1(t) 1 < (1  - ^)J  f (t-u)dG(u) 


Hence,  repeating  the  argument  of  (2.27)  - (2.31), 


(2.32)  |l3(t)-I2(t)|  < (1  - ~)2  J f(t-u)dG(2)(u), 


where 


(2.33) 


0 < a is  a constant. 


An  induction  yields  that 


CO 

(2.34)  jln(t)-I0(t)'j  < £ (1  - ®)*f*G(£)(t)  s k(t). 


0 


where  denotes  the  usual  convolution  integral. 

The  right  side  of  (2.34)  is  now  broken  up  into  a number  of  parts,  and 
upper  bounds  for  each  part  is  obtained. 


t 03 


(2.35) 


k(t)  = £ + £ (1  - 


X=0 


~)  (t) 


JL-t 


Since 


(2.36) 


if  |*G^  (t)  < K 


the  second  term  of  (2.35)  is  dominated  by 


(2-37) 


K £ (1  - f)X  < 

? t — or 

X=t 


■art 
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The  first  term  on  the  right  of  (2.35)  is  now  written  as 


(2.38)  E (1  - f)Xf*G(A}(t) 
i=0  C 

fc2  ^ a r pt-t/ in  t „t 

= 2 (1  - f)  FJ  +J  f(t-u)dGa)(u)  j. 

i=0  11  LJ0  J t-t/in  t J 


The  first  term  on  the  right  side  of  (2.38)  is  bounded  above  by 


(2.39) 


t /"(in  t) 
a V 2 


= o(t  1), 


by  the  asymptotic  behavior  of  f and  summing  the  geometric  series. 

The  second  term  on  the  right  side  of  (2.38)  is  split  into  the 
three  terms,  ignoring  distinctions  between  [c]  and  c. 


ft/ in  t t in  t t „t 

(2.40)  E + E + E J 

i=0  i=/t  /in  t i=t  in  t t-t/in  t 


f (t-u)dG^^  (u)  . 


Since  f < 1,  Chebyshev's  inequality  yields 


(2.41) 


j | f (t-u)  ]dG^  (u)  < 2 , 

t-t/in  t (t-t/in  t - ip) 


ai 


where 


(2.42) 


0 < a2  = f (t-p) 2dG(u) 
J0 


An  application  of  (2.41)  to  the  first  term  of  (2.40)  yields 


(2.43) 


f t/in  t „t 

1 £ 

i=o 


t-t/in  t 


f (t-u)dG^^  (u)  | < 


t(in  t) 


for  some  p > 0. 
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The  Central  Limit  theorem  may  be  applied  to  the  second  and  third  terms 
of  (2.40)  since  the  summation  index  i is  large. 

For  i > y t/ in  t. 


(2.44) 


G(^(t)  - 


G^Ct-t/jJn  t) 


o(- 


^)  - §(— 


a / 1 


■tf  in  t - jg|r 


a y i 


) 


Applying  (2.44)  to  the  second  term  of  (2.40)  and  using  the  mean  value 
theorem  and 


(2-45) 


9(§)  = 


yields 


(2.46) 


t In  t 

E 

f(t-u)dG(£) (u) 

M/’t/fn  t 

t-t/in  t 

< C ( t . 

-XVt/to  t9  a/£Ha(to  tVP 


< Ce-t/2  X„  t 2t!(l 
- va  in  t'5 


where  C is  a positive  constant. 

The  Central  Limit  theorem  is  applied  to  the  third  term  of  (2.40),  and 
since  the  arguments  are  large,  use  will  be  made  of  the  standard  approximation 


(2.47) 


1 - §(x)  ~ | 9(x) 


for  some  p > 0,  as  x -»  05 . 
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Applying  (2.44)  and  (2.47)  to  the  third  term  of  (2.40)  yields 


(2.48) 


t t 

S J ■ 

X=t  Xn  t t-t/Xn  t 


f(t-u)dG^  (u) 


y 


Xp.  t-Xp. 


2 ,-uVt  Xn  t 


pt  cp( 


)(at)(t/Xn  t)  -6t  Xn  t 


< 


4 2 

t p, 


< 


^e 


Xn  t 


Now  (2.37),  (2.39),  (2.43),  (2.46),  (2.48)  applied  to  (2.34)  yield 
that  for  all  sufficiently  large  t,  equation  (2.14)  holds. 


Theorem  1.  Under  assumptions  (1.1)  _to  (1.5) 


(1.6) 


lim  tP[Z(t)  > 0]  = b. 

t -+  00 


Proof . Combine  (2.12)  - (2.15)  to  yield 


(2.49) 


|p[Z(t)  >0]  - ^|  < G(n)(t)  + oCt"1). 


Let  n The  weak  law  of  large  numbers  yields  the  result. 


IIjl . Extension.  The  result  of  Theorem  1 can  be  strengthened  by  the 
method . 


Theorem  2.  Under  the  assumptions  (1.2)  - (1.5)  and  in  addition,  for  t -*  <», 


(3.1) 


t3 (l-G(t) ) - 0 


11 


and 


(3.2) 


h(3)(l)  < CO, 


then  for  t large , 

(3.3)  P[Z(t)  > 0]  ~ | -f  —|— 

t 

for  some  unspecified  constant  c. 


Remark.  Conditions  (3.1)  and  (3.2)  are  not  as  stringent  as  in  [2],  where 
more  terms  in  the  asymptotic  expansion  of  P(t)  are  given. 


Outline  of  Proof. 
Define 


(3.4) 


Jn+1 


(t) 


h(J  (t-u)dG(u) 
J0  n 


J0(t) 


t < T 

b c An  t 


t > T 


for  T » b. 

As  in  the  proof  of  (2.13)  one.  obtains 

(3.5)  0 < Un(t)  - Jn(t)  < G(n)(t). 

An  expansion  of  h about  1 to  four  terms  in  the  Taylor  expansion 
yields  that 

Jx(t)  = JQ(t)  + o(t"3). 


(3.6) 


12 


From  this  and  a similar  tedious  sequence  of  estimations  as  in  Theorem  1, 
one  obtains 


(3.7) 


|Jn(t)-J0(t)|  < G(n)(t)  + o(t“2) 


and  again  one  lets  n -»  and  applies  the  weak  law  of  large  numbers  to  complete 
the  argument. 
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